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Abstract — We provide a geometric solution to the problem of 
optimal relay positioning to maximize the multicast rate for low- 
SNR networks. The network we consider consists of a single 
source, multiple receivers and the only intermediate and locatable 
node as the relay. We construct network the hypergraph of the 
system nodes from the underlying information theoretic model 
of low-SNR regime that operates using superposition coding and 
FDMA in conjunction (which we call the "achievable hypergraph 
model"). We make the following contributions. 

1) We show that the problem of optimal relay positioning 
maximizing the multicast rate can be completely decoupled 
from the flow optimization by noticing and exploiting 
geometric properties of multicast flow. 

2) AH the flow maximizing the multicast rate is sent over at 
most two paths, in succession. The relay position depends 
on only one path (out of the two), irrespective of the number 
of receiver nodes in the system. Subsequently, we propose 
simple and efficient geometric algorithms to compute the 
optimal relay position. 

3) Finally, we show that in our model at the optimal relay 
position, the difference between the maximized multicast 
rate and the cut-set bound is minimum. 

We solve the problem for all (Ps,Pr) pairs of source and relay 
transmit powers and the path loss exponent a > 2. 

Index Terms — Low-SNR, broadcast relay channel, geometry. 

I. INTRODUCTION 

We primarily consider the problem of optimal relay posi- 
tioning in order to maximize the multicast rate in low-SNR 
networks consisting of a single source s, a set of multiple 
receivers T and an arbitrarily locatable relay r, on a 2- 
D Euclidean plane. In HI, the authors previously addressed 
this problem under a heavy and complex network flow opti- 
mization framework. They showed that optimizing the relay 
position can lead to a strong gain in the multicast rate. 

In ||2| the authors introduced equivalent hypergraph models 
for the low-SNR Broadcast (BC) and Multiple Access channels 
(MAC). The authors then derived an achievable hypergraph 
model for the broadcast relay channel (BRC), obtained by 
concatenating the equivalent BC and MAC hypergraphs. This 
concatenated model follows from constraining the source and 
relay to transmit using the optimal schemes for the low-SNR 
BC and MAC: superposition coding and frequency division, 
respectively. In this paper, building on this model, we solve 
geometrically the problem of optimal relay positioning under 



the pretext of multicast rate maximization, which is much 
simpler and efficient than the solution proposed in JT]. 

Most importantly, we establish the fact that for a given 
low-SNR BRC hypergraph Q{Af,A), the multicast rate is 
maximized by sending all the flow through at most two paths 
in succession, independently of the number of destination 
nodes. This is a consequence of simply maximizing the 
multicast min-cut. The dependency of the multicast min-cut on 
the relay position is essentially through a single path (out of 
the two), and this motivates a simple geometric interpretation 
and formulation of the problem. It should be noted that, the 
"optimal relay position" refers to the position that maximizes 
the multicast rate over a given achievable hypergraph, but in 
general the achievable hypergraph model is not necessarily 
optimal in terms of meeting the cut-set bound for low-SNR 
networks. On the other hand, the achievable hypergraph model 
performs closely to the peaky binning scheme in the case 
of a single destination fJI, and enjoys an important practi- 
cal advantage of being easily scalable to more complicated 
topologies. Finally, under our model the difference between the 
maximum multicast rate and the cut-set bound is minimized 
at the optimal relay position. 

In the proposed geometric approach, we decouple the prob- 
lem of rate maximization from the problem of computing the 
optimal relay position. This substantially reduces the complex- 
ity (compared to the flow optimization based framework in 
JT]) and also provides a great deal of insight in understanding 
the nature of such network planning problems. Finally, we 
show that at the optimal position the difference between the 
maximum multicast rate and the cut-set bound is minimized 
under the achievable hypergraph model. 

The paper is organized as follows. We introduce the low- 
SNR achievable hypergraph model of the BRC in section HI] 
Then we prove certain geometric properties of multicast in 
section [Till The computation of optimal relay position is 
divided in two parts, section |IV] for = Pr and section |V] 
for Ps ^ Pr- Finally, we conclude in section [VT] 

II. LOW-SNR SYSTEM AND HYPERGRAPH MODEL 

A. System model and notations 

The network topology is given by a hypergraph ^(A^, ^), 
where M = {s, r, T}, and all nodes except r are fixed on the 2- 
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Fig. 1. Wideband Multiple User Channels. The BC rates: Ri = (1 - p)h\ 
The BRC rates: rn = , ri 
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. Here, h gives the path loss and Dij the distance from i to j. 



D Euclidean plane. T = {ti, .., t„} denotes the set of n = |r| 
receivers ordered in increasing distance from s. C represents 
the convex hull of {s,T}. The multicast rate from s to T is 
defined as R^t = niin(i?ct), where i?»t is the total rate from 

s to receiver t £ T. Ps and Pr = jPs are the total transmit 
powers of s and r, respectively, and 7 > is their ratio. Z)„„ 
denotes the Euclidean distance between nodes u and v, and 
a > 2 the path loss exponent. For a subset Q C Af\r, define 
Lq{C) as the point in C, that minimizes the maximum over 
the distances between itself and each node in Q, i.e. 

(A) 



Ln(C) = arffmin max (Dri 

The value of objective function of the output of Program (A) 
is denoted as Dq. 

B. Low-SNR BC, MAC and BRC hypergraph models 

In H], II2J, it was shown that concatenating the low-SNR 
BC (superposition coding) and MAC (FDMA) equivalent 
hypergraph models results in an achievable hypergraph model 
for the low-SNR BRC. The rate region of this model is 
included in the capacity region of the low-SNR broadcast 
relay channel. In fact, even though superposition coding and 
FDMA are independently capacity achieving for the low-SNR 
AWGN BC and MAC channels respectively, their combination 
in general is not capacity achieving for the low-SNR relay 
channel, and a fortiori for the low-SNR BRC 

In this section, we briefly recall the equivalent hypergraph 
models for the low-SNR BC and MAC, and the achievable 
hypergraph model for the BRC [1 1. Note that in the low-SNR 
regime, BC and MAC are not limited by interference. 

1) Low-SNR BC equivalent hypergraph: Superposition 
coding is known to achieve the capacity region of the AWGN 
BC. In the low-SNR regime, the rates achieved by superpo- 
sition coding boil down to the time-sharing region ||4]-|l6l. 
For a given topology with \T\ = n receivers, the hypergraph 
will contain at most n hyperarcs with non-zero capacities IT]. 
Figures |l(a)| and |l(b)| illustrate the two-destination case. 

2) Low-SNR MAC equivalent hypergraph: In the low-SNR 
regime, interference becomes negligible with respect to the 
noise ([T], ||2l, and all sources can achieve their point-to- 
point capacity to the common destination, like with frequency 
division multiple access (FDMA). In the general wideband 
MAC with n sources, the hypergraph model consists of n 
hyperarcs of size 1 from each source Si, i £ {!,.., n} to 



the destination with non-zero capacity. Figures |l(c)| and |l(d)| 
illustrate the two-source case. 

3) Low-SNR BRC achievable hypergraph: We can obtain 
an achievable hypergraph model of the low-SNR BRC by sim- 
ply concatenating the BC and MAC equivalent hypergraphs, 
as shown in Figures [1 (e) | and [Tffl] for the two-destination case. 
As mentioned before, this achievable hypergraph model is 
suboptimal in general for the BRC, but the ability to scale 
easily to larger and complex networks is one of its biggest 
strength. 

III. GEOMETRIC PROPERTIES OF MULTICAST 
In this section, we derive the geometric properties of the 
optimal relay position maximizing the multicast rate for the 
BRC. We first focus on the single destination case of the BRC: 
the relay channel, in Section lTlI-AI Then, these preliminary ob- 
servations and properties are extended for the general problem 
with an arbitrary number of destinations, in Section IIII-BI 

A. Single destination: low-SNR relay channel 

Consider the simple network in Figure |2] (a), with a fixed 
source s, a fixed receiver t and an arbitrarily positionable relay 
r, where the multicast rate Rst from s to t is to be maximized. 
Naturally, Rst depends on the position of r. The achievable 
hypergraph in Figure |2] (a) can be broken into two subgraphs, 
shown in Figures |2] (b) and (c), which are essentially the two 
disjoint paths from s to t. 

Our claim is that the optimal position of the relay maximiz- 
ing the multicast rate from s to t lies on the line segment s — t 
joining s and t, and at this optimal position all the flow Rst is 
sent through a single path consisting of two hyperarcs, namely 
{(s,r), (r, t)} shown in Figure |2] (c). This holds true for any 
given pair of power constraints (P,, Pr) y and for any path 
loss exponent a > 2. We prove this claim in Lemmas [T] and 
ID hereafter 

We first recall the following lemma from lH]. 

Lemma 1 (Lemma 1 [T^): The optimal position of r maxi- 
mizing RsT lies inside the convex hull C. 

Here, Lemma [T] simply implies that the optimal position of 
r lies on the segment s — t. 

The rates over the three hyperarcs {(s,r), {r,t), {s,rt)} = 
A are given by. 
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Fig. 2. (a): One receiver case decomposed into two subgraphs from s to t, 
(b) and (c), respectively, (d): Optimal position of r for Ps = P,. and o = 2, 
which is at the perpendicular bisector (red) of line segment s ~ t. (e): Left 
bias for Ps < Pr- (f): Right bias for Pa > Pr- 

where A^o is the noise power spectral density. Note that the 
multicast rate is given by Rst = Rsrt + min(i?sr, Rrt)- 

Lemma 2: The optimal location of r on the segment s — t 
for a simple BRC with 7 e (0, 00) and a > 2 that maximizes 
the multicast rate Rst satisfies, 
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and the optimal (maximized) multicast rate is given by, 

Ps iPs 



(3) 



(4) 



where all the flow i?*^ is sent over the path {(s,r), (r, <)}. 

In Lemma|2]the starred entities refer the optimal values and 
for the proof the reader is referred to Appendix A. 

Lemma |2] essentially gives the position of r in terms of 
how far it is from s and r on the segment s — t. Also, it 
provides the maximized multicast rate that is achieved at 
this position. It can be easily seen that the relay position only 
affects the rate over the path {(s, r), (r, t)}. Since the min-cut 
of the path {(s, r), (r, i)} is strictly larger than the min-cut 
of the path {{s,rt)}, i.e. the rate that can be sent for a unit 
power over the former path is strictly larger than the latter path 
(Rsrt < min(i?sr, Rrt)), the rate over the path {(s, r), (r, t)} 
should be maximized first by simply maximizing its min- 
cut nim{Rsr , Rrt) before allocating any power to the path 
{{s,rt)}. The min-cut mm{Rsr,Rrt) is maximized at the 
position on the segment s — t such that rates over the two 
hyperarcs of the path {{s,r), {r,t)} become equal, and all 
the flow from s to t is transmitted over this path only. The 
maximized multicast flow is then simply given by the 
rates of either of the two hyperarcs. 

Several important conclusions can be drawn from Lemma |2] 
The multicast flow optimization can be separated from the 
determination of the optimal relay position that maximizes the 
multicast flow. Even if the aim is not to maximize the multicast 
flow (for instance by simply choosing not to use all the source 
and relay powers). Lemma |2] still gives the most suitable relay 
position for any feasible multicast rate Rst < R^f At the same 
time, the algorithmic style intuitive proof arguments in the 
previous paragraph indicate that upon computing the optimal 
relay position, the multicast rate maximization problem could 
be casted as a straightforward linear program resulting in a 



simple power allocation scheme maximizing the multicast rate. 
This fact will prove handy for the general case with arbitrary 
number of destinations. On the other hand, we observe the 
dependency of the optimal relay position on the constants a 
and 7. If 7 = 1 i.e. Pg — Pr, the optimal relay position is 
always at the mid-point of the segment s — t for any value of 
a > 2. When 7 7^ 1, there will be a natural bias on the optimal 
position of r either towards s or t, depending on the value of 
7. This bias will also depend on the value of a. Figure |2te) 
and 12 f) show the bias effect. 

B. Multiple destinations 

In this subsection, we extend the simple geometric insights 
developed in Section IIII-AI for a single destination to the 
general case of an arbitrary number of destinations |T| = n. 

Let us first note the following. For a given hypergraph 
Q{J\f,A), and a fixed position of r, we have at most + + 
(71) hyperarcs in the system, i.e. \A\ = 2n + 1. The former 
{n + 1) are source hyperarcs, emanating from s to the nodes 
in Af\s and the latter n are the relay hyperarcs, emanating 
from r to all T. Also, for any given position of r there always 
exist at least two paths that will span all the receiver set T, 
namely {(s,T)} (or {s,h..tn}) and {(s, Ti), (r, Ta)} (where 
r e Ti and Ti U T2 = {r, T}). 

Now, consider that each hyperarc (i, J) G ^ is associated 
with a continuous function fij{P^,DJj) : — > R, that 
is a monotonically increasing in the transmit node's power Pi 
and monotonically decreasing in the distance Dij, where Dij 
is the Euclidean distance between the transmit node i and the 
farthest receiver node j ^ J (from i) spanned by the hyperarc. 
Then the following theorem holds true. 

Theorem 1: Given a hypergraph Q{J\f,A) and the associ- 



ated rate functions fij{P^,D, 
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R for each hy- 



perarc in A, at the optimal position maximizing the multicast 
rate Rst one of the two multicast flow characteristics holds: 

(i) all the optimal flow i?*^ goes through at most two paths 
{(s,ri, (r,T2)} and {{s,T)}, in succession. 

(ii) all the optimal flow R*rp can be arbitrarily split between 
the two paths {(s,T)} and {(s,Ti), {r,T2)). 

For the proof of Theorem [T] refer to Appendix B. 

Theorem [U partially generalizes Lemma |2] We say partially, 
because on one hand. Theorem [T| establishes the important 
multicast flow characteristics at the optimal relay position, 
but it does not provide a simple numerical result that de- 
termines the optimal relay location (like Lemma |2). Note 
that, for a given relay position there could be multiple paths 
from s, through r, to all T, but in the Theorem [T| by path 
{(sjTi), {r,T2)} we mean the path from s, through r, to all 
T that has the highest min-cut among all the paths from s, 
through r, to all T. Intuitively, Theorem[T] states that only those 
paths will contain the multicast flow from s to the receiver set 
T that serve afl T, namely {(s, T)} and {(s, Ti), (r, T2)}. All 
other path that serve proper subsets of T will carry no flow as 
they do not contribute to the multicast flow and among all the 
paths serving all T through r, only the path with the highest 



min-cut will carry the multicast flow. This fact is a simple yet 
fundamental consequence of the definition of multicast. 

Theorem [T] reveals a lot about the nature of multicast flow 
over a hypergraph. The dependence of relay position on the 
rate of only a single path {{s, Ti), {r, T2)} reduces the problem 
to its core by removing the clutter away. In other words, now 
we only need to worry about the maximization of the flow 
over this single path and the relay position that maximizes 
the flow over this path also maximizes the multicast flow 
RsT- This result of Theorem [T| motivates a pure geometric 
interpretation of the problem. If we imagine the two hyperarcs 
(s,Ti) and (j^, T2) to be two circles Cg and Cr centered at s 
and r with radii tt^ and tt^, respectively, then the optimal relay 
positioning problem could be stated as: For a given Gi.^f,A), 
find the point in C such that when r is positioned at this point, 
max( ^77rs,7rr) is minimized while r € Cs and the region of 
union of two circles Cy = Cg U Cr encompasses all T. 

At first, it seems plausible to try a simple (preferably 
convex) optimization framework to compute such a point, 
but the condition that the two circles must encompass all M 
brings in discreteness, which we avoid for obvious reasons. In 
contrast, we propose a simple (polynomial time) algorithm to 
compute such point in the next sections. Once the optimal relay 
position is obtained, obtaining optimal power allocations (for 
s and r) maximizing the multicast rate boils down to solving 
a simple linear program involving only two paths. We divide 
the development of this algorithm into two cases of 7 = 1 and 
7 € (0, 00). The case of 7 = 1 is easy to understand and holds 
importance in its own right. In addition it develops the basic 
intuition for the proposed algorithm and leaves the extension 
to the case of all values of 7 e (0, 00), as straightforward. 

IV. [Ps = Pr) - CASE AND ALGORITHM 

In this section, we have 7 = 1 and a > 2 for a 
given GiM^A) on the 2-D Euclidean plane. The optimal 
relay positioning problem stated geometrically in the previous 
section simply boils down to finding the point in C such that 
max(7rs,7rr) is minimized while r e Cs and Cy encompasses 
aU T. We divide the problem in the following two cases based 
on the topology of the given t/(A/^, ^). 

A. s — tn mid-point case 

Lemma 3: If r is placed at the mid-point of s — tn such 
that the hyperarcs and each with radii span all 

T, then it is the optimal relay position maximizing Rst- 

The proof of Lemma |3] is a straightforward generalization 
of Lemma |2] and therefore is omitted. Intuitively, Lemma |3] 
simply states that since the farthest node (from s) t„ is also 
the limiting node for maximizing Rst, if the rate is maximized 
only to tn while guaranteeing it to all other nodes in T, then 
this maximizes Rst as well. This means that if r is placed at 
the mid-point of the segment s—tn (as this position maximizes 
the rate to tn only) and if the two hyperarcs of the path 
{(s,r), (r, i„)} {{Cs,Cr}) span all T, then clearly this is the 
relay position that maximizes Rst- 



B. General Case 

In this case we tackle all topologies and case A becomes a 
special case of it. Recall that, the entity Lq{C) represents the 
coordinates of the point which is the argument of the objective 
function of the output of program (A), and Dq is the value 
of the objective function of the output of program (A). 

Optimal relay positioning Algorithm (ORP) 

Given: g{Af,A). 

1) Compute Iq — i{A^\r}(C) and build the set No ~ {t E 
T\Dst < Di^t^Di^-^t > Dsio} = {ti,..,t^} in increasing 
order of distance from s. If No = {0}, declare Iq as the 
optimal relay position and quit, else go to step 2. 

2) Build the set Ni = {7V\(r, No)} and compute the 
point li = Lni(C). Form the hyperarcs Cg and C;^ 
of radii Dsi^ and Z?Ni, respectively. If Cy = Cg U C;^ 
encompasses all T, output li as the optimal relay position 
and quit, else go to step 3. 

3) Reform the hyperarc C^ of radius Dgt' and build the set 
N2 = {t G T\Dst > Dst'^} and compute h = LnsCC). 
Declare I2 as the optimal relay position and quit. 

Algorithm ORP is a straightforward set of basic and intuitive 
computational steps based on the properties of the point 
^0 = Lj^\r{C). If there exist no node t' G T such that 
t' ^ Cs and Dst' < Di^t' (i-e. set Nq is empty), that can be 
directly reached by s rather than by a path through r, then is 
certainly the optimal relay position. In contrast, if the set Nq 
is not empty, then there exist at least one receiver node in the 
system that influences the computation of the optimal relay 
position but can be served directly by Cg. Therefore, either 
the nodes in No can be removed from the computation of 
the optimal relay position {li in Step 2) and max(7rs, tt^) can 
be further reduced or we could reform the hyperarc Cs with 
radius i'st;^ (where, t\n is the farthest node in Nq from s) and 
then computing the point I2 for the nodes that were not covered 
by Cs and thus reducing the value of max(7rs, tt^). Note that. 
Algorithm ORP categorizes all possible topologies of the given 
g [Af, A) in three steps and there is no underlying iterative 
process. This makes algorithm ORP behave like a numerical 
formula, which we originally wanted from Theorem [T] 

We leave the formal proof that ORP always outputs the 
optimal relay position maximizing Rst to Appendix C and 
extend this simple approach in a straightforward manner to 
the case of all values of 7 G (0, 00) in the next section. 

V. Ps ^ Pr- CASE AND ALGORITHM 

In this section, we consider 7 £ (0, 00) for a given 
g{Af,A) and a > 2. Almost all the theory developed in 
Section |IV] simply transcends to this section, with certain 
notable differences. Mainly, that when 7 7^ 1 it gives rise 
to a bias in the positioning of r ( ref. Figure |2e) and |2jf)). 
Taking into account the bias while computing the optimal 
relay position will be the main enhancement in this section. 
Likewise previously, we first consider the s — tn case. 



A. s — tn case 

Lemma 4: Given Q{M ^A), if r is placed on s — t„ at a 
distance of Dsr = x^^^ from s, such that r £ Cg and Cy = 
Cs U Cr spans all T, then it is optimal relay position that 
maximizes Rst- 

The line of argument for the proof of Lemma |3] (using 
Lemma |2|i could be simply generalized for Lemma |4] 

B. General Case 

In this case, like in Section IIVI we generalize to all 
topologies. As we know, that the values of 7 (when not equal 
to 1) and a inflict the bias on the relay position. The main 
difference in case of Pg ^ Pr is the computation of the point 
k = LQ{C) (i = {0,l}), given by, 

k = Lq{C) = argmin ( max ( ^Dg^, Dij) I . (B) 
«ec \U&Q\s) ) 

and the computation of the set Nq = {i G Tj ^Dsia > 
Di„f} = {t[,..,t'^}, in the Algorithm GRP. Program (B) 
and the set Nq takes into account the bias induced by the 
differences in the transmit power of the source and relay and 
the value of a. The rest of the algorithm remains the same. 

Now that we have an efficient algorithm for computing the 
optimal relay position, we can be more ambitious to assess 
the standing of our work in a more theoretical sense. Gne 
of the important consequences of this work that signifies its 
theoretical importance is shown in Figure |3] We computed the 
difference between the optimal multicast rate R*rp (for a given 
position of r) and the cut set bound for |r| = 9 receiver nodes 
network at 21 interesting positions, including the optimal relay 
position computed by the Algorithm GRP. At the optimal relay 
position (blue point), this difference is minimized, confirming 
the fact that the optimal relay position not only results in gains 
but the maximized multicast rate is theoretically closest to the 
cut-set bound at the optimal relay position in our framework. 

It is worth mentioning that the theory developed in this 
paper well transcends to the low-SNR fading channels , which 
we do not discuss here but can be easily generaUzed from the 
results of E) and ID. 

VI. CGNCLUSIGN 

We list the important deductions from our work in the 
following points. 

1) The problem of optimal relay positioning to maximize 
the multicast rate for the achievable hypergraph model 
of low-SNR networks using superposition coding and 
FDMA, can be decoupled from flow optimization and 
casted as a simpler geometric problem, as opposed to a 
complex network optimization approach of [1]. 

2) The geometric properties of multicast are innately simple 
and provide interesting insights for relay positioning 
problem. This is largely due to the fact that all the 
multicast flow is pushed over at most two paths which 
is a direct consequence of the definition of the multicast 
flow, and this results in simple geometric interpretation. 




Fig. 3. |T| = 9 case with green receivers, red source and blue as tlie optimal 
relay position. The optimal -R^t and cut set bound dilference (in nats/sec) 
is calculated for 21 positions and is the lowest at the optimal relay position 
(blue). We assume = = 1 (normalized) and a = 4. 

3) Importantly, the benefits of determining the optimal relay 
position are substantiated by the fact that the difference 
between the maximized multicast rate and the cut-set 
bound at the optimal position is minimized. 
We now outline, what we think are certain important future 
directions our work could take. The geometric properties of 
multicast give great insights and are surprisingly easy to work 
with. This motivates us to ask further, whether is it possible 
to apply the simple techniques of our work for the optimal 
relay positioning problem to moderate and high-SNR regimes 
that are interference limited. Another natural and interesting 
dimension is to look at the possibility of extending this work 
to multicommodity flows. 
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Appendix A 
Proof of Lemma |2] 

Proof: We only consider the positions in the interior of 
the segment s — t. Then, the multicast rate is given by 

Rst = Rsrt + min(i?sr7 Rrt) 
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By assumption, we have D^t > max(Dsr, £'rt)- Thus, in 
the minimization of Q, the first and the second term are 
respectively decreasing and increasing affine functions of A. 
Two cases can occur If ^Dgr > Drt, then the second term 
is always larger than the first term, which consequently is the 
minimum of the two. The first term decreases in A, thus Rst is 
maximized for A = 0. Else, if ^Dgr < D^t, the two affine 
functions intersect in the interval [0, 1] at A = 1 — The 
multicast rate Rst is maximized at this intersection. Note that 
for the position of r satisfying ^Dsr — Drt, both solutions 
match: A = 1 - = 0. 

By Lemma [U the relay position maximizing the multicast 
rate lies on segment s — t. Then, we can write 



Ds 



Ds 



Dr 



(6) 



and the relay position is simply determined by the distance 
Dsr- Using (|6]l, the conditions "^i^Dsr ^ £'rt in Lemma |2] 
can be rewritten in function of Z),,. as 
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Ds 



Given the optimal power allocation A*, and using (|7]), the 
multicast rate Rst can be rewritten as the following function 

of D,r 



R. 
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if n < _Dst_. 
if n > -Pst 



(8) 

From ds), it can be seen that Rst is an increasing function of 
Dsr over (0, j^^], and then a decreasing function of Dsr 
over [j^^^, Dst)- Therefore, the multicast rate is maximized 
at the border between these two intervals: D* = , 'f^'V - 
Substituting D*^ in dSJl yields i?*^. This completes the proof 
of the Lemma. ■ 

Appendix B 
Proof of Theorem |7] 

A simple assimilation of the basic graph theoretic and 
Euclidean geometric concepts helps form the fundamental 
reasoning for the proof. Let's assume that the hypergraph 
tJ(A/', .4) is given with the constant 7 G (0,oo) (where, 
R,. = 7P5) and each hyperarc {i, J) G A is associated with 
any continuous rate function fij{P^ ,D~j) : — > IR, 
that is monotonically increasing in the transmit power of 



the emanating node i of the hyperarc and is monotonically 
decreasing in the distance Du (between the transmit node i 
and the farthest node j G J from i). We notice that there 
are only two transmitters in the system s and r and the 
multicast rate from s to the receiver set T is defined as 
RsT = min(i?5t), where Rst is the total rate received by 

the receiver t € T. 

Proof: To prove the theorem, we first notice that for a 
given position of r there are at least two paths, namely certain 
paths of the type {(s, Ti), {r, T2)} (where, Ti U T2 = T) and 
the path {(s,T)}, that span the whole receiver set T. Any 
other path, that only serves the proper subsets of T does not 
count in contribution to the multicast rate Rst- 

Among all the paths from s to T, that go through r (i.e. of 
the type {{s,Ti), (r,T2)}, where, T1UT2 = T), only the path 
with highest min-cut contributes to the multicast flow Rst- 
Let us denote this path as {(s, T{), ir, T^)}, where T[\JT!^ = 
T. Once the min-cut of the path {(s, T{), (r, Tj)} is reached, 
considering it has the highest min-cut among the paths that 
span all T through r, no flow can be sent over any other path 
of the type {(s, Ti), (r, T2)}. This is true because when the 
min-cut of the path {(s,T{), {r^T!^)} is achieved (for a fixed 
position of r) either Ps is consumed or Pr is consumed. If 
Ps is consumed, no more multicast flow can be pushed, and 
if Pr is consumed before Ps then rest of the flow have to 
be pushed over the path {(.s,T)} (not involving r). On the 
other hand for a given position of r, if the min-cut of the path 
{(s,r{), (r,r^)} is strictly less than of {(s,T)}, flien afl the 
multicast flow must be sent over the path {(s,T)}. 

This implies that for any given position of r, all the 
multicast flow must be sent over at most these two paths. 
Now, we can write down the min-cut of the multicast flow as 

RsT ^fsT{P+,D;^) + 

mm{fsTi{p+,D;^,), frTi{P+, d;^,)). 

Now, consider the region Cn = Cs D Ct^, which is the 
intersection of the two circles centered at s and r with radii 



min{Ds 



2D, 



and TTr ~ min(£'s 



2 <^Ds 



respectively. The radii tt^ and tt^ takes the bias due to a and 
7 into account. Simply stated, if 7 > 1 then tTs < TTr, and if 
7 < 1 then TTs > TTr, and finally if 7 = 1 then the two circles 
have equal radii. It is clear that if 7 e (0, 00) then the area of 
Cn > 0. 

If the relay is positioned outside Cn, then 
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2^Dr 



1+ ^' 1 



> Ds 



implying, 

fsT{P^,D:^) > min(/,T,(P+,I):T,)./rT.(P+, A"tJ) • 

This means that the min-cut of the path {(s,T)} is strictly 
larger than the min-cut of the path {(s, T[), {r, T2)}, implying 
that all the multicast flow must be sent over the path {(s, T)}, 
rendering relay useless. Hence, the optimal relay position must 
lie inside Cp. 
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Fig. 4. T = {t-i,t2,tz} case illustrating step 3 of the algorithm ORP. (a): 
The relay r is placed at the point Zq = Lj^\^^{C) with No = fsti < 

Drtj and Dsti > Dsr (b): Reforming the hyperarc Cs and r is placed at 
I2 = (C) (where, N2 = {t2,t^}), thus reducing max(7rs, TTr). 

From here, it is straightforward to see that if the optimal 
relay position lies in the interior of Cn, rendering the min- 
cut of the path {(s,T{), (r, T2)} strictly larger than the path 
{(s,T)}; then flow over the path {{s,T[),{r,Tl^)} must be 
maximized first and then the flow over the path {(s,r)}, in 
order to maximize Rst- This proves the first point of the 
theorem. 

Similarly, if the optimal relay position lies on the boundary 
of the region Cn, then 

fsT (P+ , D;^) = min {fsn (P+ , ^.T, ) , frT. (Pr : ^."t. )) : 

rendering the min-cut of the two paths equal. In this case, all 
the flow can be sent over the path {(s, T)}, {(s, Ti), (r, T2)} 
by arbitrarily sharing the flow between them. This case is 
reminiscent to the case when relay is placed outside Cn, 
but for completeness we count it as an individual case, and 
moreover in this case relay is not really useless. This proves 
the second part and hence completes the proof of Theorem |3] 



Appendix C 
Proof of optimality of Algorithm ORP 

Assume a given G{J^,A) and 7 = 1. The argument of the 
output of Program (A) is a point I and the objective function 
value of the output of Program (A) is distance denoted by Dq 
(where, Q is the set of points of input to Program (A)). 

Proof: In order to prove that Algorithm ORP always 
outputs optimal relay position, we need to prove that the 
three steps suffice to tackle all the topologies of a given 
Q{J\f,A) (namely, the distribution of the points Af\r on the 
2-D Euclidean plane). 

First, we divide all the topologies in two classes. In the first, 
the point / is the optimal relay position (which corresponds to 
the step 1 of the algorithm ORP), and the second class in 
which the point / is not the optimal relay position (this class 
corresponds to the Step 2 and 3 of the algorithm ORP). The 
only complicated case (if at all) is the Step 3, so we will go 
about proving the optimality of the output of algorithm ORP 
backwards in the order Step 3, then Step 2 and finally Step 1. 



For a given Q{J\f,A), compute the point Iq defined as, 

lo ~ SiTg min (max(_Dij)) 

subject to: i £ C. 

Form the hyperarcs Cs of radius Dsi„ and Cr of radius 
Dj^\^. Denote the value of the quantity max(7rs,7rr) ~ (■ 
Construct the set No = {t G T\Dst < Di„t & Di„t > 
Dsig} — {t'l, in the increasing order of distance from 

s. Considering the set No is not empty, take the farthest node 
tj„ from s. If Dst'^ > Dgig, then it is clear that t'^ should 
be approached directly from s and not through r, because 
> Dsi„ and ma.x{Dsio, Di^t'^) < Dst'^. Reforming a 
source hyperarc Cs of radius Dst'^ (where, Dst'^ < C)- the 
set N2 = {t G T\Dst > Dst' } can be constructed consisting 
of all the nodes not lying in the area Cs- Now, computing 
the point I2 = LnsCC) we could form the second hyperarc 
Cr of radius -Dn2- Note that -Dn2 < C because the set N2 
consists only the nodes in T that are not in the hyperarc Cs 
and Dst'^ > Dsi„- Denoting niax(7r5, tt^) = (" (with respect 
to point I2), we now have C," < (. We cannot further reduce 
max{TTs,TTr) as tj„ is the farthest node in No that satisfies 
this property. Thus I2 is the optimal relay position. Figure |4] 
illustrates this step for |r| = 3 case. 

On the other hand, if the node tj„ satisfies the relation 
Dst' < Dsio^ it is clear that all the nodes in No could 
be dropped from the computation of the point I and the set 
Ni ~ {J\f\{r, Nq)} can be constructed. Therefore, computing 
the point li = Lni(C), gives the optimal relay position as 
there is no node in the set Ni that influences the computation 
of the optimal relay position unnecessarily. Again, reforming 
the hyperarcs and denoting max(7rs,7rr) ~ C (with respect 
to point li), we can easily see that (' < (. The value of 
max(7rs,7rr) cannot be reduced further because there is no 
receiver node in T that is in the set No that cannot be 
encompassed by the area of union of the two hyperarcs Cs 
and Cr (constructed with respect to the point ^i). Thus, in 
this case li is the optimal relay position. 

Finally, if the set Nq = {0}, the point Iq is clearly the 
optimal relay position as there is no receiver node in the 
system that is affecting the computation of the relay position 
and can be dropped off simultaneously. Hence, the three steps 
of algorithm ORP always outputs the optimal relay position 
for a given Q{JV,A). 

■ 

The case when 7 7^ 1 is a straightforward generalization 
and line of argument for the proof of optimality remains the 
same for the case of 7 = 1. 



